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1. Introduction 

The various evidences of the existence of dark matter have led to the development of exper- 
iments dedicated to detect dark matter directly. The detection relies on the measurements of the 
recoil of atoms hit by a dark matter candidate. One popular class of dark matter models involve an 
interaction between a WIMP and a Nucleon mediated by a Higgs exchange. Therefore, the scalar 
content of the nucleon is a fundamental ingredient in the WIMP-Nucleon cross section. In this 
way, the uncertainties of the scalar content translates directly into the accuracy of the constraints 
on beyond the standard model physics. Since the coupling of the Higgs to quarks is proportional to 
the quark masses, it is important to know how large scalar matrix elements of the nucleon are, in 
particular for the strange and charm quarks. 

One common way to write the parameters entering the relevant cross section are the so-called 
sigma- terms of the nucleon: 

OnN = ni{N\uu + dd\N) and OQ = m{N\uu + dd-lss\N), (1.1) 

where m denotes the light quark mass. Quantifying the scalar strangeness content of the nucleon a 
parameter yj^ is introduced, 

_ 2{N\ss\N) 
{N\uu + dd\N) 

which can be also related to the sigma terms of the nucleon in eq. ( |l.l[ ). 

The direct computation of the above matrix elements is known to be challenging on the lattice 
for several reasons. First, it involves the computation of "singlet" or "disconnected" diagrams 
that are very noisy on the lattice. Second, discretisations that break chiral symmetry generally 
suffer from a mixing under renormalization between the light and strange sector, which is difficult 
to treat in a fully non-perturbative way. However, twisted mass fermions offer two advantages 
here: they provide both an efficient variance noise reduction for disconnected diagrams and a 
convenient way to avoid the chirally violating contributions that are responsible for the mixing 
under renormalization. 



2. Lattice Setup 

In our simulations we use the mass-degenerate twisted mass action in the light sector and 
the mass non-degenerate twisted mass action in the strange and charm sector. The quark masses 
of the heavy quark doublet have been tuned such that the Kaon and D-mesons masses assume 
approximately their physical value. The reader interested in more details about aspects of this 
setup is referred to ||l], §]. The twisted mass action in the light sector reads 
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where the hopping parameter K = (2amo + 8r)~ is defined in terms of amo, the bare Wilson mass, 
r is the Wilson parameter and /x^ is the bare twisted mass parameter. The Wilson parameter is fixed 
to |r| = 1 is all our simulations. When K is tuned to its critical value a situation called maximal 
twist is achieved which guarantees 0{a) improvement of physical observables. 

For further needs we also introduce the operators ± denoting the upper and lower compo- 
nents of the Wilson twisted mass operator in flavour space (also referred to as the Osterwalder- 
Seiler Dirac operator): 

D,,±[U]=tr^^D^[U], (2.2) 

where tr denotes the trace in flavour space. Dq^±[U] then corresponds to 1-flavour twisted mass 
operators with Wilson parameter r = ±1. 

When we discuss the 2-point and 3-point functions necessary for this work, we will use the 
so-called physical basis of quark fields denoted as Xj/q. This field basis is related to the twisted 
quark field basis, Xq by the following field rotation 

xifq = e''^r^''x, and xif^ = x/'^r^'\ (2.3) 

where the twist angle Wi = n/lat maximal twist. 

In order to compute matrix elements involving strange quarks, we will work within a mixed 
action setup. For reasons that will become clear later, we choose to introduce in the valence sector 
an additional doublet of degenerate twisted mass quark of mass ^iq . The mass fiq can then be 
tuned to reproduce the Kaon and D-mesons mass in the unitary setup. Preliminary estimates of the 
matching masses gives ajJ-g = 0.0185 in the strange sector and ajj.,. = 0.2514 in the charm sector 
and we will approximately use these values for [iq further on. In this contribution, we work at 
a fixed lattice spacing corresponding to a 0.078 fm with mpsL > 4 and pion masses ranging 
approximately from 300 to 500 MeV. 

3. Matrix elements 

In the following, Yq with index l,s,d will denote the quark fields of the light (l), strange (s) or 
charm (c) quarks in the physical basis. In order to be self-contained, we describe in this section the 
relevant correlation functions used in this work. The nucleon two-point function reads : 

Cw,2pts(^ -^src^src) = I^tr r±(^^(x)^(Xsrc)), (3.1) 

X 

where Xsrc = (^ro^src) is the source point and the subscript refers to the proton or to the neutron 
states for which the interpolating fields are given by: 

= e"'" (if^'^'^Ysd^) u' and = e"^' [d'^^'^'^Ysu^) d' . 

The projectors used are = and is the charge conjugation matrix. 
The nucleon three-point functions is 

4?pts(f^'^^op,Xsrc) = L tr r±(^A'WO,(Xop)7^(Xsrc)), (3.2) 
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where Oq is an operator having scalar quantum numbers, e.g. Oq = qq, Afop = fop — hvc is the time 
of insertion of the operator, and = t — gives the so-called source-sink separation. Note that in 
the twisted basis the scalar operators read 

Oq = iXqYi'C^Xq, whcrc q = l,s,c , (3.3) 

and are hence given by the pseudo scalar density. The flavour structure of the operators in the 
twisted basis will be crucial in the following. 

Since we consider an operator with a non- vanishing vacuum expectation value, we also define 

^A'jSpts (^ii'^^opi-^src) = 3p''jj,(f5, Afop,.Xsrc) — 2pts(^)-^src) ^(0(y(Xop)) . (3.4) 

The desired scalar matrix elements can then be computed using the asymptotic behaviour of the 
ratio of a three and two-point functions: 

Ro,{t.,Atop) = = i^M^) + ^i^'^'") + ^{e-^^'^-'"^^) . (3.5) 

^A',2ptsV^-^srcj 



The general form of the 3-point functions in Eq. (3.2) lead to both, connected (C, illustrated in 



fig. [T]a) and disconnected illustrated in [T]b) contributions, 

^iV,3p'ts(^*''^^op) ~ ^A',3pte(^-* ''^^op) "I" ' (^.V) ^opj-^src) (3.6) 

In the following we will denote by /?conn. (resp. /?disc.) the contribution of C^'^' (resp. to the 



ratio defined in Eq. (3.5). The sum of the connected and disconnected contribution to the ratio will 
be denoted Rfa\\. 

In order to improve the signal over noise ratio, we have averaged the disconnected part over 
forward and backward propagating proton and neutron states. In addition, we have used up to 4 
randomly chosen source points per configuration for the 2-point function computation to enhance 
our statistics. 

As will be detailed in a forthcoming publication |^], the operators Oq do not mix under renor- 
malization and hence have a straightforward renormalization pattern very similar to chirally invari- 
ant overlap fermions. We consider this fact as a major advantage of our twisted mass approach for 
computing the scalar quark contents of the nucleon. 

I 
I 
I 

t 

I 

I t 





(a) (b) 
Figure 1: We illustrate the connected (left) and the disconnected (right) graphs that arise from the contrac- 
tions leading to the 3-point function discussed in the text. 
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4. Variance reduction 



The main challenge to compute the 3-point functions of Eq. ( [3.2[ ) is to evaluate . Our 

strategy is based on a variance reduction technique for twisted mass fermions introduced in [Q] 
and used to study the r\' meson in ||^. It relies on the fact that in the twisted basis the disconnected 
contribution that has to be evaluated is related to the difference of 1 /Dq,^ — 1 

Here, we will employ the one-end-trick [Q] to compute the disconnected distribution stochas- 
tically evaluating 



lia^qY. 0[r]W75r(^>[,](x) =£trr 



1 



1 



[X,X) + I 



(4.1) 



where 



(/)[,] = (l/D,,+ )^[,] and 0[;] = ^[;](1/D,,_ 



(4.2) 



where we have introduce Nr independent random volume sources denoted For the generation 
of the random sources we have used a Z2 noise setting all field components randomly from the set 
{1,-1}- 

In our tests for the signal to noise ratio (SNR) we first investigated how the SNR depends on the 
number of stochastic sources Nr used. We found that for Nr>7 there is no significant improvement 
of the SNR. Nevertheless, we have used 12 stochastic sources per configurations in all our results. 
In Fig. ^ we compare the SNR of the twisted mass specific variance reduction technique to a more 
standard method based on the the hopping parameter expansion ||5|]. We show the ratio Rq^ of 
eq. (3^), for a fixed value of A?op = t^/l = 6 as a function of the number of configurations N. We 
conclude that the SNR is increased by a factor ~ 3 with our improved noise reduction technique 
which allows to obtain a result at the 5 a significance level with only a moderate statistics. 





Hopping parameter expansion 
• twisted mass noise reduction 






5o limit 


V* - ♦ ♦ ♦ 

* ♦* ♦ ♦ 




♦ ♦♦ • 








I0 limit 




100 200 300 400 


500 600 


700 



• Flo< 
A Rdi 

• R|y 



A ^ ^ A 



15 



Figure 2: Signal to noise ration (SNR) of the ratio 
Ro,j for a fixed time Atop = ts/^ = 6 as a function of 
the number of gauge field configurations A^. Rp^ is 
evaluated here in the strange quark regime. We com- 
pare our noise reduction technique, with the hopping 
parameter expansion technique. The dashed lines in- 
dicate the 1(7 and 5(7 significance levels and the short 
dotted lines are only shown to guide the eyes. 



Figure 3: Plot of the two contributions Rdisc. and 
^conn. to the ratio Rq^ relevant for the extraction of 
(^kn- Rm\ is the sum of both contributions. 
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5. Light (7— term 



In Fig. ^ we show the results obtained for the bare ratio Rf^w introduced in Eq. (3^) from 
which OnN can be computed. The connected part, /?conn.» for a source-sink separation of = 
12a, is shown by the black filled circles. It has been computed using the standard method of 
"sequential inversions through the sink". As can be seen, 7?conn. shows a time dependence indicating 
excited state contributions, a systematic effect that has not been taken into account yet in this work. 
The disconnected part, 7?disc., is represented by blue triangles in Fig. ^ The disconnected part is 
significantly smaller than the connected part /?conn. and contributes at the ^ 10% level to the full 
ratio /?fuii represented by the red diamonds. 

We have computed o^n so far for three values of the pion mass, see table Our present data 
do not allow for a reliable extrapolation to the physical point for which additional data at more pion 
masses would be necessary. 



mps (MeV) 


OnN (MeV) 


318 


99(6) 


392 


152(9) 


455 


228(15) 



Table 1: Fit results for o-ji^ as a function of the pion mass. Only statistical error are estimated. 



6. Strangeness of the nucleon 

In Fig. we show /?disc. for a quark mass of ajXq = 0.018 corresponding approximately to 
the strange quark mass. The in principle freely selectable source-sink separation has been fixed 



to 12 lattice units. The ratio Rq^ of eq. ( |3.5[ ) shows a time dependence indicating that also in the 



case of the strange quark excited states may be important. We nevertheless extract a plateau value 
as indicated in Fig. ^ which is clearly different from zero. Combining this value with the result 
for the scalar matrix element obtained in the light quark sector discussed above, allows us finally 
to compute y^. We have performed such a computation at four values of the pion mass at fixed 
value of the lattice spacing and the results are summarized in Fig. ^. Performing a simple linear 
extrapolation we obtain, as shown by a red star, = 0.069(27) in the chiral limit, where only 
the statistical errors have been taken into account. Although an investigation of systematic effects 
are still missing, we checked that by varying the bare strange quark mass to apLg = 0.016 does not 
change significantly the value of yj^i . 

As a final remark we mention that we have also computed the charm quark content of the 
nucleon. Unfortunately, here the SNR is of order one and hence no clear signal can be extracted. 
From our present data we can only provide a qualitative estimate that {N\Oc\N) < {N\Os\N). 

7. Conclusion 

In this proceedings contribution we have shown that with twisted mass fermions it is possible 
address the disconnected contribution to the scalar matrix elements of the nucleon. This becomes 
especially important when the strange and the charm content are computed since there only discon- 
nected graphs appear. In addition, with twisted mass fermions the renormalization pattern of such 
matrix elements is the same as for chirally invariant discretizations. As a result we were able to 
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Figure 4: Time dependence of /?disc. in the strange 
quark mass regime (a/X^ = 0.018). The source- 
sink separation has been fixed to 12a (842 config- 
urations) 



Figure 5: Our data foryAf as a function of the pion 
mass. The data are extrapolated to the chiral limit 
using a simple Unear extrapolation. 



obtain accurate results in the light and strange sector at fixed lattice spacing and for several quark 
masses. Our main result is a value yf^ = 0.069(27). This value is compatible with recent lattice 
results obtained by several groups [|8|, ^ |l^. The still missing systematic uncertainties on this result 
will be addressed in future simulations. 
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